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Abstract 

We measure by Monte Carlo simulations ^/string for a model of random surfaces 
embedded in three dimensional Euclidean space-time. The action of the string is 
the usual Polyakov action plus an extrinsic curvature term. The system undergoes 
a phase transition at a finite value A c of the extrinsic curvature coupling and at the 
transition point the numerically measured value of 7 s t r m 3 (A c ) 0.27 ± 0.06. This is 
consistent with 7 s t r m 3 (A c ) = 1/4, i.e. equal to the first of the non-trivial values of 
^/string between and 1/2. 
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Introduction 



Since the first studies of regularized string theory, i.e. the theory of random surfaces, 
it has been clear that bounds on the string susceptibility exponent 7 (also called 
the entropy exponent) made it difficult to obtain interesting models, i.e. models 
which have interesting continuum limits when the cut-off is removed. In almost all 
models considered so far 7 < 1/2 is a rigorous bound [1, 2]. In the simplest random 
surface models like the hyper-cubic lattice model where one sums over all planar 
random surfaces built of plaquettes of the lattice, it is in addition possible to show 
that 7 > implies 7 = 1/2. The value 7 = 1/2 is the value for surfaces which are 
degenerated to so-called branched polymers, and this limit is uninteresting from a 
string theoretical point of view. If we consider string theories as special examples of 
matter fields coupled to 2d quantum gravity we have in this broader context many 
examples of exactly solvable models where 7 can be calculated analytically. The 
fractal structure of the random surfaces is a function of the central charge of the 
conformal field theories coupled to 2d quantum gravity [3] : 



However, as is seen from (1) 7 < in all cases where the formula makes sense, 
i.e. for c < 1. So far there have been no exactly solvable models with c > 1 and 
this includes unfortunately the bosonic string theories of interest since they have 
c = d } the dimension of target space. It has been difficult to obtain genuine random 
surface models or 2d gravity models coupled to unitary matter theories 5 with 7 > 0. 
Recently it has been understood that if more than one coupling constant is present, 
as will be the case for instance if we consider multi-Ising models coupled to gravity, 
it is possible to obtain a non-trivial behaviour where 7 = l/(n + l),n>2 [7] (see 
also [8]). However, the corresponding behaviour will be related to the behaviour of 
a spin system with central charge c < 1 (i.e c = 1 — and 7 = — 1/n) and 



The simplest case is c = and leads to 7(c) = 1/3. This is realized in the case of 
infinitely many Ising models coupled to gravity [9, 10], and can be considered as a 

5 It should be noted that it for some time has been known how to penetrate the barrier 7 = 
by modifying the concept of random surfaces, such that not only piece-wise linear manifolds but 
more general simplicial complexes were allowed. In this way one could obtain 7 = 1/n by using 
generalizations of the multicritical matrix models, but the corresponding c would be less that zero, 
i.e. non-unitary [4, 5, 6]. 
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mean field result. However, for a specific multi-spin model it is not known which c 
is selected and this is still a major gap in our understanding of the back-reaction 
of matter on 2d quantum gravity. Nevertheless, numerical simulations indicate that 
the back-reaction for relatively small c (c < 4) is a universal function of c, like for 
c < 1, i.e. the fractal geometry of the random surfaces is the same whether we 
couple eight Ising spins or we couple four q=4 states Potts models to 2d quantum 
gravity [11]. 

If we turn from 2d quantum gravity to the genuine random surface models we 
can enlarge the coupling constant space by adding an extrinsic curvature term. The 
motivation for this is that it can be shown that the string tension does not scale in 
the models without such a term. This is true both for the hyper-cubic model [1] 
and the model defined by summing over triangulated surfaces [12]. In both cases it 
is believed that the non-scaling of the string tension is a consequence of the domi- 
nance of branched polymers which again might be related to the tachyon problem 
of the bosonic string. The dominance of branched polymers will diminish if we add 
an extrinsic curvature term. Further, it can be argued that one gets an extrinsic 
curvature term by integrating out the fermions of the superstring [13]. For these 
reasons the models of random surfaces with extrinsic curvature have been exten- 
sively studied over the last couple of years in the context of dynamical triangulated 
random surface models. The present letter reports on the measurement of 7 at the 
critical point where the bare string tension and the bare mass in the model scales 
to zero. 



The model 

The model of random surfaces with extrinsic curvature we will discuss has the fol- 
lowing partition function [14, 15, 16, 17]: 

N 

Z(A, A) = £ / n<**? * (i) (£4) e" 5T , (3) 

TET N J i=l 

where T denotes a (closed) triangulation of spherical topology, consisting of of N 
vertices, T/v denotes the class of all such triangulations, xf the space-time coordinate 
associated with vertex i of the triangulation, and // the space-time index (// = 
1,2,3 for d = 3). Finally St denotes the action corresponding to the particular 
triangulation: 

St = ^EK "4) 2 + A (l" cos %)). (4) 
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In this formula (ij) is the link connecting the vertices i and j (if any) and 6^ 
the angle between the normals to the two oriented triangles (ijk) and (jil) with a 
common link (ij). This action is the discretized version of the following continuum 
action: 

S(\) = J d 2 T^jg af3 (d a x^dpx^ + \d a n^dpn^) + m 2 J <P Ty /g (5) 

The first term is the bosonic string term in the formulation of Polyakov and the 
A-term is a scale invariant extrinsic curvature term. In (5) g a p refers to the intrinsic 
metric. A related version of this action is 

S(\) = J d 2 rVhh af3 ( i i d a x^d p x^ + \d a n^dpn^), (6) 

where h a p refers to the induced metric: 

dx^ dx^ 

Since d a n^ = H a ph f3l d 1 x fl (Weingarten eq. for the normals in terms of the second 
fundamental form H a p) eq. (6) can be written in a way often used in theory of fluid 
membranes: 

S(\) = J dA(p + H 2 ), H = Q- + 1) , H 2 = H aP H lS h^h rjS (8) 

where H denotes the extrinsic curvature and r\(r) and r 2 (r) are the principle radii 
of curvature at the point x a1 (ti,t 2 ). (5) and (6) are not entirely equivalent at the 
classical level, but they both produce smoother surfaces with increasing A. The 
action in the form (7) or (8) is not easy to use in the discretized approach since care 
has to be exercised in order that each physically distinct surface is counted only 
once. This is automatically satisfied if we use (3) and (4). 

Let us briefly summarize what is known about the phase structure of the model. 
At A = we have a regularized version of Polyakov's bosonic string theory in d = 3. 
This theory is known to have a non- vanishing bare string tension [12]. For increasing 
A we get smoother strings at short distances, but the string still seems to belong 
to the same universality class as the pure bosonic string. For A — > A c 1.45 we 
have a genuine phase transition [15, 16, 17, 18] to a string theory where both the 
mass gap and the string tension seem to scale [17, 19]. It is at this point where we 
measure the susceptibility exponent 7(A C ) which is also called the entropy exponent 
since it is related to the finite volume partition function. In the discretized approach 
the internal area of the surface (the equivalent of the continuum / d 2 T \/g) is simply 
proportional to the number of triangles Nt } since each triangle is assumed equilateral 
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in the internal metric. We will take the lattice units such that the area is just 

J d 2 r^ = A N = N T = 2N - 4 (9) 

where N is the number of vertices of the closed surface of spherical topology. For a 
fixed number of vertices, or equivalent a fixed internal area A, we have 

Z(X, A) ~ A^~ 3 e^ A (l + 0(-^j . (10) 

This formula provides us with a convenient way of measuring 7(A) by counting so- 
called baby universes. A minimal neck baby universe, abbreviated "minbu" is a part 
of the surface connected to the rest of the surface by a minimal loop of links, i.e. a 
loop of length three. From (9) it follows that the average number of such minbus of 
area a < A/2 for surfaces with the fixed total area A is 

N A (a) ~ (a(A-a)yW- 2 . (11) 

It is important to note that the leading, non-universal part of the distribution 
Z(X } A) is not present in Na(cl), and this makes it a convenient observable to use for 
the determination of 7(A). This was first pointed out in [20] and has been verified 
in practice in a number papers, even in higher-dimensional gravitational theories 
[21, 22, 11]. 



Numerical results 

In the simulations of the model we used the vectorization over number of systems 
proposed in [25] and then developed in [26]. Most of the measurements were per- 
formed on lattices of the sizes 800 and 3200 triangles. They were oriented on collect- 
ing histograms for the baby universe distribution. To estimate the autocorrelation 
time we additionally measured the radius of gyration which is known to be the 
slowest mode in the algorithm. On the smaller lattice we simulated 64 systems in 
parallel, while for the larger, because of the memory limitations we had to lower 
their number to 32. The measurements were taken every 100 sweeps. The measure- 
ment of the minbu distribution is not vectorizable, so we ran it in the scalar mode. 
The maximal size of the lattice in our simulation was not only limited by the mem- 
ory but also by the autocorrelation time r which is particularly high in the critical 
region where we spent most of the computer time. From the measurements of the 
gyration radius we estimated r. In the critical region r is 0.13(4) • 10 6 sweeps for the 
3200 triangle lattice and 0.14(3) • 10 5 sweeps for the 800 triangle one. For the larger 
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lattice we run 3.0 • 10 6 sweeps in each of 32 systems while for the smaller one 1.0 • 10 6 
sweeps in each of 64 systems. This corresponds around to 23, 70 autocorrelation 
times, respectively. At the beginning we made 10 6 and 10 5 sweeps, respectively, 
for thermalization, which in both cases corresponds to 7r. One can worry about 
insufficient thermalization. We checked, however, that discarding the data from 
the next few autocorrelation times does not change the the resulting histograms for 
minbu distribution. In fact, the autocorrelation times for the quantities describing 
the intrinsic geometry are much lower for the studied range of lattice sizes than for 
the gyration radius [23]. 

We simulated systems for different values of the coupling A. The resulting minbu 
distribution can be divided into three classes corresponding to two different phases 
of the model and to the critical region. In fig. 1 we show typical curves for each of 
these cases for the system with 3200 triangles. They are plotted in log log scale. 
Going downwards the curves correspond to A = 0.5, 1.45, 1.9. 

The slopes of the curves in the crumpled phase and in the critical region differ 
slightly. Both are close to —2. As one can see from the formula (11) the slopes 
are dominated by the gamma independent power —2 in the minbu distribution. To 
visualize the 7-dependent information in the minbu distribution, we plot in fig. 2 
the function F(a) = A^(a)(a(l — a/A)) 2 which explicitly removes the bulk power 
—2. The function F(a) is plotted in log log scale against ail — a/ A). It is clearly seen 
that the upper curve, A = 0.5 has a negative slope, although it's value changes along 
the curve. This is a well known problem for 2d gravity models coupled to purely 
gaussian fields or multi-Ising spins with c > 1 [11] and might be due to logarithmic 
corrections. We obtained this type of curves for small A's which correspond to being 
deep in the crumpled phase, and it confirms that the model is dominated by the 
gaussian term and belongs to the same universality class for this range of A. The 
gaussian model itself was discussed in more detail elsewhere [11]. 

When approaching the critical region the curves gradually change until in the 
critical region their slope becomes positive. The situation is in this respect very 
similar to that observed in the spin models where the effective numerical value of 
7, as a result of the finite size effects, passes gradually from that for pure gravity 
outside the critical region, to the value characteristic for critical matter sector, in the 
vicinity of the pseudocritical value of the coupling constant. In our case the critical 
region, where the slope stayed roughly constant corresponds to A between 1.4 and 
1.5. As values of A for which we performed the high statistics simulation we chose 
A = 1.44 for 800 triangles and A = 1.45 for 3200. The choice of the pseudocritical 
values of the coupling constant was dictated by estimation of positions of the heat 
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N t 


7 


800 
3200 


0.31(7) 
0.27(6) 



Table 1: 

capacity maxima from the results for the smaller systems and is slightly higher than 
the value A c 1.425 quoted in [18]. As was explained above the small difference 
has no effect on the measured value of 7. 

The lowest curve in fig. 1 corresponds to the 'flat' phase where the normals become 
almost aligned. In this phase the acceptance rate goes down below 10% and the 
autocorrelation time r becomes even larger. The phase seems to be a lattice artifact 
and one would need much larger lattices and probably a new algorithm in this phase. 

To extract 7 we fitted our data points directly to the formula (11). As input to 
our fits we took the data points Na(cl) averaged over all systems. Their errors SNa(cl) 
were naturally provided by the widths of the distributions in different systems. As is 
clearly seen from the figure 2, both a quality of the fit and the value of 7 depend on 
the range of the minbu areas a used for the fit. We made different fits for different 
ranges a mt -„ < a < a max . The values of 7 are plotted for different values of the 
lower limit a mt -„ in fig. 3 and fig. 4. Different symbols correspond to different values 
of the upper limit a max . The fitted values of 7 saturate with increasing a mt -„. The 
corresponding values of % 2 decrease. At a mt -„ around 30, x 2 becomes of the order one 
per degree of freedom and from there on it stays at this level only weakly depending 
on a max . This means that when one excludes the minbus of the size smaller than 
30 one gets surprisingly good fit to the formula Na(cl) (11) in the broad range of a 
without including the next-to-leading terms which were necessary, for example, in 
the pure gravity model. Although the values of 7 look very stable as a function of 
the limits a mt -„ and a max one can observe a tendency that they seem to decrease when 
the limits are increased. This is seen in fig. 4 for the 3200 triangle lattice, where the 
effective 7 goes gradually down after the maximum around a mt -„ = 30. The values 
decreased also with increasing a max . The same tendency is also reflected in the fact 
that the values of 7 are lower for the larger lattice. All this can be summarized by 
saying that the effective value of 7 decreases for larger areas. The results for the 
effective 7 are summarized in the table 1. 
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Discussion 

The measurements performed in this paper were done for systems with finite size. 
The measured values of 7 C indicate that 1/4 < 7 C < 1/3. The error bars quoted 
are probably over-estimated. For increasing size we see a clear trend to lower the 
value of 7 C which suggests that for the infinite system we have 7(A) 1/4 for the 
ensemble of random surfaces described by the partition function (3) at A A c . 
As explained in the introduction this implies that the surface theory given by (3) 
realizes the first non-trivial exponent 7 = 1/4 of the possible values 7 = 1/n, 
n > 2. The corresponding "effective" central charge is c = 1/2 and it indicates an 
equivalent fermionic description of the random surfaces with extrinsic curvature. It 
is worth to recall that the same observation has been made a long time ago for hyper- 
cubic surfaces in d=4: Extensive computer simulations [24] showed that 7 1/4 
for a model where "back-tracking" of plaquettes 6 was not allowed. On the other 
hand 7 1/2 for d=8 in the hyper-cubic lattice model without back-tracking. This 
suggests that at least for low dimensions we have a universal non-trivial susceptibility 
exponent for random surfaces with extrinsic curvature. It should be most interesting 
to check at which dimension the exponent 7 of dynamically triangulated surfaces 
changes from 1/4 to 1/2, and to get an understanding of this from the point of view 
of continuum string or field theory. 
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Figure captions 



Fig.l The minbu distribution Na(ci) in log log scale, for the lattice with A = 3200 
triangles, for the different values of the coupling A : 0.5 (+), 1.45 (O), and 1.9 

(°) 

Fig. 2 The function F(A) = Na(ci) ■ (a(a — 1/A)) 2 versus a in log log scale for the 
lattice with A = 3200 triangles. 

Fig. 3 The effective values of 7 for the lattice with A = 800 triangles, for different fit 
ranges : a mt -„ runs along the vertical axis, while different a max are represented 
by : 200 (O), 400 (+), 600 (□) and 800 (x). 

Fig. 4 The effective values of 7 for the lattice with A = 3200 triangles, for different fit 
ranges : a mt -„ runs along the vertical axis, while different a max are represented 
by : 150 (O), 200 (+), 250 (□) and 300 (x). 
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